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Abstract 


Severe acute respiratory syndrome (SARS) is a respiratory disease caused by a newly found virus, called SARS coronavirus. In 
this study, the cleavage mechanism of the SARS coronavirus main proteinase (M?*° or 3CLP*°) on the octapeptide NH>- 
AVLQ J SGFR-COOH was investigated using molecular mechanics and quantum mechanics simulations based on the experimental 
structure of the proteinase. It has been observed that the catalytic dyad (His-41/Cys-145) site between domains I and IJ attracts the x 
electron density from the peptide bond Gln-Ser, increasing the positive charge on C(CO) of Gln and the negative charge on N(NH) 
of Ser, so as to weaken the Gln—Ser peptide bond. The catalytic functional group is the imidazole group of His-41 and the S in Cys- 
145. Ng on the imidazole ring plays the acid—base catalytic role. Based on the “distorted key theory” [K.C. Chou, Anal. Biochem. 
233 (1996) 1-14], the possibility to convert the octapeptide to a competent inhibitor has been studied. It has been found that the 
chemical bond between Gln and Ser will become much stronger and no longer cleavable by the SARS enzyme after either changing 
the carbonyl group CO of Gln to CH, or CF, or changing the NH of Ser to CH, or CF3. The octapeptide thus modified might 
become an effective inhibitor or a potential drug candidate against SARS. 
© 2004 Elsevier Inc. All rights reserved. 
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A new coronavirus (SARS CoV)! has been found as 
the etiological agent of the disease SARS [1-4]. SARS 
CoV is a positive-sense, single-stranded RNA virus with 
the largest viral RNA genome known to date [5,6]. The 
SARS CoV genome is composed of around 29,700 
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nucleotides. The replicase gene alone encompasses more 
than 21,000 nucleotides, 2/3 of the whole genome se- 
quence. The replicase gene encodes two overlapping 
polyproteins, ppla (4377 a.a.) and pplab (7072 a.a.). 
SARS CoV main proteinase (M°"°) cleaves the polypro- 
tein at no less than 11 conserved sites, a process initiated 
by the enzyme’s own autolytic cleavage from ppla and 
pplab [7,8]. The released polypeptides mediate all of 
the functions required for SARS viral replication and 
transcription. The functional importance of the MP*° 
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in the viral life cycle makes it an attractive target for the 
drugs design against SARS [7-10]. 

Anand et al. [7] suggested that the rhinovirus 3C?"° 
inhibitor AG7088 could well serve as a starting point 
for an anti-SARS drug based on the theoretical homol- 
ogy model of SARS CoV M?"®. Chou et al. [9] found the 
fitting problem of AG7088 to the binding pocket of 
SARS CoV M"°, and they suggested its derivative 
KZ7088 as a better starting point. Sirois et al. [11] did 
virtual screening for SARS-CoV protease based on 
KZ7088 pharmacophore points. Chou et al. [9] also sug- 
gested an octapeptide NH>-AVLQ J SGFR-COOH that 
can be converted to an inhibitor of SARS CoV MP’® by 
some chemical modification. This octapeptide was syn- 
thesized, proved to be a cleavable peptide, and showed 
strong inhibiting activity in a tube test [12]. Du et al. 
[13] analyzed the cleavage mechanism of the octapeptide 
by SARS CoV M?"° and discussed the possible chemical 
modification of the octapeptide using molecular dynam- 
ical and quantum mechanical simulations based on the 
homology structure of SARS CoV M?"®, to convert it 
into a stable and effective drug. This progress on finding 
inhibitors against the SARS enzyme is also summarized 
in a recent review [10]. 

The experimental structure of SARS CoV M?"° was 
first determined by Rao’s research group [8]. They pro- 
vided several crystal structures of SARS CoV M?"° at 
different pH conditions and the complex structure of 
SARS CoV M?*° with ligand hexapeptidyl CMK [7,8]. 
The experimental structure and the theoretical homol- 
ogy structure of SARS CoV MP” share many structural 
similarities, such as the catalytic active region and the 
Cys-His cleavage dyad [7,8]. However, there are several 
differences between the theoretical and the experimental 
structures. For examples, the experimental structure of 
SARS CoV MP"® is a dimer, the protomer A has cata- 
lytic activity, and the protomer B has no activity. How- 
ever, the N finger of protomer B, which plugs in the 
active cleft of protomer A, plays an important role in 
both dimerization and maintenance of the active form 
of the enzyme [8]. These new findings from the experi- 
mental structure of SARS CoV MP’ provide a reliable 
basis for inhibitor design. 


Method 


The octapeptide NH»-AVLQ J SGFR-COOH shows 
strong inhibiting activity in vitro [12]. However, the use 
of native peptide for clinical applications has been ham- 
pered by their intrinsic metabolic instability (indigestion 
by enzymes) and poor absorption through membrane. A 
cleavable peptide of SARS CoV M?"° means that it pos- 
sesses a competent binding with its receptor and a 
breakable peptide bond in the cleavage site. To convert 
a Cleavable peptide into a stable and effective inhibitor, 


the peptide has to be modified according to the “‘distort- 
ed key”’ theory [14,15]. One important step for the chem- 
ical modification is a good understanding of the 
cleavage mechanism, particularly the changes in chemi- 
cal bonds [16-22]. 

The protease-susceptible sites in a given protein or 
peptide usually extend to an octapeptide region [14,15]. 
The corresponding amino acid residues are sequentially 
symbolized by eight subsites Ry, R3, Ro, Ry, Ry, Ry, 
Ry, Ry and the eight matching positions in the protease 
by Sy, S3, S2, Sy, Sy, Sy, S3, Sy, respectively (see Fig. 3 
of [15]). Occasionally, the susceptible sites in some pro- 
teins may contain one subsite less or one more. How- 
ever, eight amino acid residues are the most common 
case. Although the protein being cleaved contains much 
more than eight amino acid residues, only an octapep- 
tide segment fits in the active region of the protease. 
Therefore, our research will focus on the cleavability 
of octapeptides. 

As shown in Fig. 1, the cleavage point is always on 
the peptide bond between R; and Ry». In Fig. 1A we 
use the combination of a thin line and a dashed line to 
represent the conjugate m property in the peptide bond. 
In Fig. 1B the breakable peptide bond between R, and 
R, is converted to a strong “hybrid peptide bond” 
through a chemical modification and hence cannot be 
cleaved by the protease. The modified octapeptide can 
still bind to the active site, but it cannot be cleaved by 
the enzyme, just like a distorted key that can neither 
open the lock nor be pulled out. The octapeptide thus 
modified will automatically become a desired inhibitor 
candidate. 

The octapeptide AVLQSGFR is the first suggested [9] 
based on the molecular structure of SARS CoV M?P"° 
and has been proved cleavable by experiments. In this 
research we study the cleavage mechanism, the proper- 
ties of the relevant chemical bonds, and the catalytic 
interactions between the octapeptides and the SARS 
CoV MP® using molecular mechanical and quantum 
chemical simulations to provide useful insights for the 
chemical modification. The study was carried out in 
the following four steps: (1) docking the octapeptide to 
SARS CoV MP"® followed by an energy minimization 
for the complex structure using molecular mechanics, 
(2) computing the atomic charge distribution in the 
binding pocket for the energy-minimized structure using 
ab initio quantum mechanics, (3) computing the molec- 
ular energy, chemical bond properties, and atomic 
charges of the octapeptide with the background charge 
distribution [23] of SARS CoV MP*® using ab initio 
quantum chemistry, and (4) with the same background 
charge distribution, computing the molecular energy, 
chemical bond properties, and atomic charges of the 
modified octapeptide. The computations have been car- 
ried out with the molecular dynamics and standard ab 
initio quantum chemistry at the HF/6-31G* level [23]. 
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A 


SARS Cov MP? 


SARS Cov MP’? 


Fig. 1. Schematic illustration to show (A) a cleavable octapeptide is chemically effectively bound to the active site of SARS enzyme and (B) although 
still bound to the active site, the peptide has lost its cleavability after its scissile bond was modified from a hybrid peptide bond to a strong bond. 


Adapted from Chou [15] with permission. 


Results 
Docking and energy minimization 


The experimental structure (lujl.pdb) of SARS CoV 
MP" [8] is used in the calculation and shown in Fig. 2A. 
Each protomer structure of the SARS CoV M?"° dimer 
has three domains. The two protomers A and B form a 
right angle and the join part is in the domain III. The N 


INA&B 


N-Finger 
of B 


inhibitor 


Protomer B Protomer A 


finger of protomer B plugs into the domain II of pro- 
tomer A and plays an important role in the catalytic 
activity. The PDB data of SARS CoV M?"° contains 
300 hydration water molecules, which provide the solva- 
tion ion environment with some water molecules joining 
the catalytic reaction. In Fig. 2A the green signs ‘‘<”’ re- 
present the solvation water molecules. The docking 
calculation between SARS CoV M?"® and the octapep- 
tide AVLQSGFR was performed using the molecular 


Octapeptide 


lil A&B 


Fig. 2. (A) Crystal structure of SARS CoV MP"® consists of two protomers A and B. Each protomer has three domains and the two protomers form 
a right angle at the domain III. The N finger of protomer B plugs into the domain I of protomer A. The green signs ““<”’ represent the solvation 
water molecules. (B) The complex structure obtained by docking the octapeptide to SARS CoV M?"° and a follow-up energy minimization. The 
catalytic dyad of His-41 and Cys-145 is marked by an arrow and the peptide bond to be cleaved of the octapeptide is shown in green. 
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operating environment program package [24]. A total of 
25 docking conformations were obtained, and the one 
with the lowest energy was used for further energy min- 
imizing calculation. Fig. 2B shows the complex structure 
of SARS CoV M?"° with the ligand octapeptide. The 
positions of catalytic dyad His-41 and Cys-145 are indi- 
cated in Fig. 2B. 

According to Anand et al. [7], Yang et al. [8], and 
Chou et al. [9], the catalytic active region is in the pocket 
between domains I and I of protomer A. During calcu- 
lation the catalytic active region on which we were 
focusing contains the amino acid residues Cys-22, Gly- 
23, Thr-24, Thr-25, Leu-27, His-41, Val-42, Cys-44, 
Thr-45, Ala-46, Glu-47, Asp-48, Met-49, Leu-50, Asn- 
51, Pro-52, Tyr-54, Phe-140, Leu-141, Asn-142, Gly- 
143, Ser-144, Cys-145, Glu-163, His-164, Met-165, 
Glu-166, His-172, Asp-187, Arg-188, and Gln-189 from 
promoter A and the residues Phe-3, Arg-4, Lys-5, Met- 
6, and Ala-7 from promoter B. But in the practical cal- 
culation, the actual region involved contains much more 
than the above-listed amino acid residues. It includes all 
the amino acid residues within 5A surrounding the 
octapeptide [9,18]. 

Fig. 3A shows the catalytic active region of SARS 
CoV M?"® and the location of octapeptide (white). The 
catalytic dyad of His-41 and Cys-145 is shown by a ball 
and stick drawing, the peptide bond Gln-Ser is in green, 
and the N finger in pink. Fig. 3B gives a close look at the 
bridge water molecule (W) between the SARS CoV MP"° 
and the octapeptide. The distance O(W)—C(CO) is 
4.70 A, H,(W)— Ns)(His-41) is 2.39 A, Ho(W)—-N(NH) 
is 4.20 A, and H(NH)-S(Cys-145) is 4.12 A. The bridge 


Gln—Ser 
bond 


Cys-145 


N-Finger of B 


water molecule plays an important role in the proteolytic 
reaction. These results are derived from the complex 
structure obtained by following the aforementioned 
docking and energy-minimization procedures. 


Hydrophobic-hydrophilic interaction 


The hydrophilic and hydrophobic complementary 
interaction is important in the ligand-receptor interac- 
tion. According to the « sphere model [25,26], the active 
region of proteinase consists of a series of « spheres. 
Each interaction site between the ligand and receptor 
contains one or more o spheres. There must be at least 
one hydrophobic interaction site. Figs. 4A and B are 
the hydrophilic and hydrophobic surfaces of octapeptide 
AVLQSGER and the catalytic active region of SARS 
CoV MP"®, respectively. In the figures the green parts 
are hydrophobic surfaces and the blue parts are hydro- 
philic surfaces. There is good hydrophilic and hydro- 
phobic complementary interaction between ligand and 
receptor: the hydrophilic parts of the ligand correspond 
to the hydrophilic parts of the receptor and the hydro- 
phobic parts of the ligand correspond to the hydro- 
phobic parts of the receptor. The octapeptide 
NH>-AVLQ J SGFR-COOH has four hydrophobic 
amino acid residues: R4(Ala), R3(Val), R3’(Phe), and 
R4(Arg). Among them R4(Ala) and R3(Val) are ex- 
posed in solvent and R3’/(Phe) and R4’(Arg) are covered 
by hydrophobic surfaces of the proteinase. The peptide 
bond to be cleaved between R1(Gln) and R1‘(Ser) is 
just in the join point between the exposed part and the 
covered part. 


Fig. 3. (A) Catalytic active region of SARS CoV M?” and the location of octapeptide (white). The catalytic dyad of His-41 and Cys-145 is shown in 
ball and stick form, the peptide bond Gln-Ser is in green, and the N finger is in pink. (B) Bridge water molecule (W) between SARS CoV MP’® and 
octapeptide. The distance O(W)—C(CO) is 4.70 A, H)(W)—Ngi(His-41) is 2.39 A, H2(W)-N(NH) is 4.20 A, and H(NH)-S(Cys-145) is 4.12 A. 
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Ala \R4) 


Phe (R37) 


Val (R3) 


Arg (R4’) 


Fig. 4. (A) Hydrophilic and hydrophobic surfaces of octapeptide NH-AVLQ | SGFR-COOH. The green parts are hydrophobic surfaces and the 
blue parts are hydrophilic surfaces. The four hydrophobic amino acid residues are R4(Ala), R3(Val), R3’(Phe), and R4’(Arg). (B) Hydrophilic and 
hydrophobic surfaces of SARS CoV M?"°. The hydrophobic residues R4(Ala) and R3(Val) are exposed in solvent, while R3’(Phe) and R4’(Arg) are 


covered by hydrophobic surfaces of the proteinase. 


Quantum chemical calculations 


The calculations were focused on the effect of SARS 
CoV M?*° on the peptide bond to be cleaved. How to 
properly treat a biomacromolecule and the solvent ef- 
fects is a nontrivial problem in quantum chemical calcu- 
lations. Electrostatic interaction plays the dominant role 
in ligand-receptor combinations and must be taken into 
account in the quantum mechanical study. In view of 
this, we treated the SARS CoV M?"® as a background 
with many point charges. The amino acid residues in 
the catalytic region between domains I and II of pro- 
tomer A were divided into seven segments, and their 
atomic charges were separately calculated using the ab 
initio quantum chemical HF/6-31G* method. The 5 
amino acid residues of the N finger of protomer B were 
taken into account also. The 76 amino acid residues in 
eight segments are listed in Table 1. There are 1148 
atomic charges in the background, including all atoms 
in the catalytic cleft; the overlapped atoms are not 
counted. 

Another difficult problem for quantum mechanical 
calculation for a biological system in solution is the sol- 
vent effect. In the PDB data [8] of SARS CoV M?"® there 


are 300 solvation water molecules. These water mole- 
cules not only provide the solvation environment but 
also indicate that some water molecules are involved 
in the catalytic reaction. In our calculation the 300 sol- 
vation water molecules were treated as the background 
charges. Each water molecule has 3 point charges. 
Therefore, there are 1148 + 300 x 3 = 2048 background 
point charges. 

In Table 2 we list the Mulliken atomic charges gM" 
and electrostatic potential equivalent charges q**? of 
His-41, which are obtained by fitting atomic charges to 
the electrostatic potential at the van der Waals surface. 
The atomic charges from semiempirical method AM1 
are quite different from the results by ab initio HF/6- 
31G* [23]; in particular; the g*°? charges of AM1 do 
not appear reasonable. Because the HF/6-31G* q®*? 
charges reproduce the quantum mechanical electrostatic 
potential on the molecular surface, in this research we 
use g®°” from the HF/6-31G* calculations to illustrate 
our points. The imidazole group of His-41 is a nucleo- 
philic agent and its pXk=6.0, near the concentration 
of proton [H“] in pure water. The imidazole ring is both 
a proton donor and a receptor. It has good catalytic 
activity in the physiological condition [27]. It can be seen 


eee of amino acid residues in the active cleft of SARS-CoV M?"®: 76 amino acid residues and 1148 atoms are included 

Domain I of A Domain II of A B 
Al-l A 1-2 A 1-3 A 2-1 A 2-2 A 2-3 A 2-4 Bl-1 
20-30 40-50 49-59 140-150 160-169 168-178 181-190 1-5 
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Table 2 
Atomic charges of amino acid residue His-41 of SARS CoV MP", in 
atomic unit e 


Atom AMI AMI HF/6-31G HF/6-31G 
qa" gh? ia gee 
N —0.5314 0.5569 —0.9054 —0.4429 
H(N) 0.2987 —0.5387 0.4113 0.1614 
Cy —0.05334 2.0820 —0.0462 0.1192 
Hy, 0.2611 —0.7777 0.3081 0.2461 
Cc 0.3174 —1.9381 0.7638 0.3254 
O —0.3257 0.0303 —0.5750 —0.5518 
Cp —0.2302 —1.2490 —0.3088 —0.1595 
Hg3 0.2633 —0.1937 0.3268 0.2057 
Hp 0.1628 0.1624 0.1903 0.0576 
C, —0.0590 1.7519 0.2424 0.2770 
Ns —0.2861 —1.7978 —0.8863 —0.4094 
Hs1(+) 0.4132 0.7184 0.4821 0.4318 
Cs —0.1669 1.5835 0.1266 0.0852 
Hys2 0.3390 —1.5772 0.3350 0.0992 
Car 0.0476 0.1000 0.5776 0.2349 
Hey 0.3489 1.6428 0.3430 0.2178 
No —0.2709 —3.4887 —0.9061 —0.4587 
Her 0.4020 2.1826 0.46957 0.4604 


from Table 2 that the negative charge of Ns; (—0.4094) 
and the positive charge of H;, (0.4318) are large in mag- 
nitude and are close to the peptide bond Gln-—Ser. 
Therefore, the proton Hj, attracts the n electron density 
from the peptide bond Gln—-Ser so as to weaken this 
chemical bond. The polar nitrogen N35, attracts one 


hydrogen H, of bridge water and makes the water mol- 
ecule dissociated, causing the nucleophilic attack of hy- 
droxyl group OH to the C(CO) and the electrophilic 
attack of H* to N(NH) on the two sides of peptide bond 
Gln-Ser. The sulfur atom S in Cys-145 forms a hydro- 
gen bond with H(NH) and helps the proteolytic reaction 
(Fig. 3B). 


Chemical modification of cleavable peptides 


The peptide bond is considered a pseudo m bond. The 
six atoms in the peptide bond (Gln)Ca CO-NHCa (Ser) 
are almost in a plane. Fig. 5A shows the electron density 
contour map of the peptide bond Gln-Ser in the gaseous 
phase on the z plane containing the carbonyl group CO 
of Gln and the NH group of Ser. The electron density 
contours surrounding atoms N(NH) and C(CO) form 
two triangles reminiscent of the hybrid orbits of the 
sp’ type, and the 2p, orbits, which are perpendicular 
to the plane and form a 7} bond system. In Table 3, 
we list the calculated atomic charges of the six atoms 
on both sides of the peptide bond Gln-Ser obtained 
from the ab initio HF/6-31G* calculations in the gas- 
eous phase with the background charges of SARS- 
CoV M""°. The negative charge of N(NH) of serine in 
the octapeptide increases to —0.8612 in the background 
of proteinase charges from the original value of —0.8534 
in the gaseous phase. Similarly, the corresponding 


T T 
12 12.5 13 13.5 14 14.5 15 15.5 16 16.5 17 17.5 


Fig. 5. (A) Electron density contour map of the peptide bond Gln—Ser in the x plane containing the carbonyl C and O atoms of Gln and N(NH) of 
Ser, as calculated in the gaseous phase. (B) Contour map of the differential electron density of peptide bond Gln—Ser in octapeptide AVLQSGFR, 
obtained by subtracting the electron density in the gaseous phase from the electron density in the proteinase background and solvent water. 
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rie charges of relevant atoms in peptide bond Gln-—Ser of octapeptide, in atomic unit e 
Phase Gln side Ser side 

Co C(CO) O(CO) Ca N(NH) H(NH) 
g’>?(Gas) —0.0133 0.8534 —0.6252 0.4124 —0.8050 0.3528 
gs? (MP"’) 0.0020 0.8612 —0.6219 0.3842 —0.8190 0.3497 


positive charge of the carbonyl carbon C(CO) in gluta- 
mine of the octapeptide increases to 0.8190 from the 
original 0.8050. This is favorable to the electrophilic 
attack of OH™ to C(CO) and the nucleophilic attack 
of H’ to N(NH) in the proteolytic reaction. 

Fig. 5B is the contour map of differential electronic 
density of the peptide bond Gln-Ser in the octapeptide 
AVLQSGER, obtained by subtracting the electron den- 
sity in the gaseous phase from the electron density in the 
background [23] of SARS CoV M?"° and solvent water 
molecules. In Fig. 5B the bold gray line is the 0 value 
line where the electron densities in the gaseous phase 
and in the protease background are equal. The solid thin 
lines show the regions where the computed electron den- 
sities are greater in the SARS CoV M?"° background 
than in the gaseous phase, and the dashed thin lines 
show the areas where computed electron densities are 
smaller in the proteinase background than in the gas- 
eous phase. We find that along the peptide bond be- 
tween (Gln)C—N(Ser) the electron densities increase on 
the N(Ser) side and decrease on the C(GlIn) side. This 
change is favorable for a nucleophilic attack of anion 
OH to (Gln)C and for an electrophilic attack of cation 
H* to N(Ser). 

SARS CoV M?"° has very high cleavage specificity 
and is highly conserved in the coronavirus family [7]. 
The cleavage specificity of SARS CoV MP"® concen- 
trates on the sites R2, Rl, and R1’. The development 
of peptides as clinically useful drugs is limited by their 
poor stability and bioavailability, which is due in part 
to their inability to readily cross membrane barriers such 
as the intestinal and blood-brain barriers. Systematic 
chemical modification strategies that convert peptides 


into drugs are an attractive research topic in current 
medicinal chemistry [28]. For the peptide inhibitor of 
proteinase, the chemical modification to cleavable octa- 
peptide should focus on the scissile peptide bond be- 
tween R, and R, to be cleaved by SARS CoV M?”. 
A successful example shows that after the peptide bond 
CO=NH is replaced by reduced single bond CH,—NH, 
the modified octapeptide rennin inhibitors have shown 
greatly increased affinity for rennin and strong resis- 
tance to enzyme hydrolysis [29]. 

In this study we suggest a similar modification to 
octapeptide AVLQSGFR. In the peptide bond (Gln) 
CO=NH(Ser) the carboxyl group CO of glutamine on 
subsite R, is replaced by CH or CF) or the NH of ser- 
ine on subsite R, is replaced by CH); then the m bond 
system is broken and the modified octapeptide AVLQ 
SGFR may become a stable inhibitor for SARS-CoV 
MF*° and an effective drug candidate against SARS. 
Here we show the possibility of the chemical modifica- 
tion through computational modeling. 

Table 4 lists the atomic charges of the six atoms on 
both sides of the hybrid peptide bond Gln-Ser after 
the chemical modification. Comparing with Table 3, 
we find that after CO is replaced by CH2, the charge 
of the carbon atom in CH, turns negative (—0.6572) 
from originally a positive value (0.8612) in the CO 
group; hence the nucleophilic attack by OH” becomes 
impossible. Furthermore, the negative charge of N(NH) 
on the Ser side decreases to —0.2892 from —0.8190 in 
Table 3, and hence the electrophilic attack by H* also 
becomes more difficult. Likewise, the modifications 
CO > CF, and NH > CH) will make the proteolytic 
reaction more difficult also. 


Table 4 
Atomic charges in peptide bond Gln-Ser after modification, in atomic unit e 

Gln side Ser side 
CO > CH, Co C(CH2) H(CH2) Co N(NH) H(NH) 
g°>?(Gas) 0.0770 —0.6756 0.1351 —0.1689 —0.2473 0.1910 
gs (MP"’) 0.0867 —0.6572 0.1251 —0.2121 —0.2892 0.1807 
CO > CF, Co C(CF3) F(CF)) Co N(NH) H(NH) 
g°>?(Gas) 0.0814 0.3679 —0.0625 0.2955 —0.6854 0.3143 
gs? (MP"’) 0.1182 0.3092 —0.0653 0.2597 —0.6483 0.3085 
NH > CH, Co C(CO) O(CO) Co C(CH2) H(CH2) 
g’>?(Gas) 0.0311 0.8581 —0.5980 0.4387 —0.6707 0.1546 
ges? (MP"’) 0.0551 0.8437 —0.5940 0.4086 —0.6460 0.1511 
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Discussion 


Searching for cleavable peptides is the first step in 
finding potential inhibitors against the proteinase. Even 
within a limited range of octapeptides, it is by no means 
an easy job because the possible different patterns of 
octapeptides are as many as 20° = 2.56 x 10'°. However, 
some algorithms, such as the vectorized sequence-cou- 
pled algorithm [14] and the discriminant function algo- 
rithm [15], can be successfully used to narrow the 
search scope. Meanwhile, similar methods that can be 
used to search the cleavable peptides by SARS enzyme 
have been developed [5,6]. 

The octapeptide NH >-Ala—Val—Leu—Gln-Ser—Gly 
Phe-Arg-COOH is the first designed on the basis of 
the molecular structure of SARS-CoV MP" [9] and it 
has been experimentally demonstrated to be a cleavable 
peptide with high bioactivity [12]. A detailed study of 
the cleavage mechanism needs a great deal of experimen- 
tal work. However, computational simulation may pro- 
vide useful insights for an in-depth understanding of the 
mechanism. The cleavable octapeptide could be changed 
into an effective inhibitor of SARS-CoV MP?" or a drug 
candidate against SARS after a proper chemical modifi- 
cation to replace the cleavable peptide bond with a non- 
cleavable bond. The octapeptide thus modified will lose 
its susceptibility although it can still tightly bind to the 
active site and become a competent inhibitor, just like 
the role of a distorted key in blocking the access of other 
keys to a lock, as elucidated by Chou [15] in a review 


paper. 
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